---
title: "P(ID) for n = `r params$n_graph`"
output: 
  flexdashboard::flex_dashboard:
    orientation: columns
    social: menu
    source_code: embed
params:
  n_graph: 10
  n_samp: 10000
---

```{r setup, include=FALSE}
library(dygraphs)
library(quantmod)
library(flexdashboard)

root <- rprojroot::find_root(".git/config")
r_dir <- file.path(root, "r")
cpp_dir <- file.path(root, "cpp")
invisible(lapply(list.files(r_dir, full.names = TRUE), source))
invisible(lapply(list.files(cpp_dir, full.names = TRUE), sourceCpp))

# grid experiment full
n_graph <- params$n_graph
res <- replicate(length(gspec), NULL)

set.seed(2023)
for (i in seq_along(gspec)) {
  
  grid <- gspec[[i]]$grid
  var_params <- names(grid)
  grid[["n"]] <- n_graph
  df <- prob_id_grid(adj_gen = gspec[[i]]$adj, cfd_gen = gspec[[i]]$cfd, 
                     grid = grid, nsamp = params$n_samp)
  res[[i]][["raw"]] <- df
  if (length(grid) == 2L) {
    
    by_args <- c("params1", "n", "adj", "cfd")
  } else by_args <- c("params1", "params2", "n", "adj", "cfd")
  res[[i]][["smry"]] <- df[, lapply(.SD, mean, na.rm = TRUE), by = by_args]
  res[[i]][["smry"]] <- res[[i]][["smry"]][complete.cases(res[[i]][["smry"]])]
  res[[i]][["vis"]] <- vis_prob_id(res[[i]][["smry"]], var_params)
}
x <- latex2exp::TeX("my title $\\alpha = \\beta$")

adj_labs <- list(Gnp = latex2exp::TeX("$G(n, p)$"),
                 unif = latex2exp::TeX("Uniform"),
                 scale_free = latex2exp::TeX("$BA(n, a, \\alpha)$"))
cfd_labs <- list(Gnq = latex2exp::TeX("$G(n, q)$"),
                        prunning = latex2exp::TeX("IP(q)"))
```

Row
-------------------------------------
    
### `r attr(adj_labs[[1]], "latex")` + `r attr(cfd_labs[[1]], "latex")`
    
```{r}
res[[1]]$vis
```
 
### `r attr(adj_labs[[1]], "latex")` + `r attr(cfd_labs[[2]], "latex")`
    
```{r}
res[[2]]$vis
``` 

Row
-------------------------------------
    
### `r attr(adj_labs[[2]], "latex")` + `r attr(cfd_labs[[1]], "latex")`
    
```{r}
res[[3]]$vis
```
    
### `r attr(adj_labs[[2]], "latex")` + `r attr(cfd_labs[[2]], "latex")`

```{r}
res[[4]]$vis
```

Row
-------------------------------------
    
### `r attr(adj_labs[[3]], "latex")` + `r attr(cfd_labs[[1]], "latex")`
    
```{r}
res[[5]]$vis
```
    
### `r attr(adj_labs[[3]], "latex")` + `r attr(cfd_labs[[2]], "latex")`

```{r}
res[[6]]$vis
```